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Dissociative ionization of : Few-cycle effect in the joint electron-ion energy spectrum 
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Joint electron-ion energy spectra for the dissociative ionization of a model HJ in few-cycle, infrared laser 
pulses are calculated via the numerical ab initio solution of the time-dependent Schrodinger equation. A strong, 
pulse-dependent modulation of the ionization probability for certain values of the protons’ kinetic energy (but al¬ 
most independent of the electron’s energy) is observed. With the help of models with frozen ions, this feature— 
which mistakenly might be attributed to vibrational excitations—is traced back to the transient population of 
electronically excited states, followed by ionization. This assertion is further corroborated employing a two- 
level model incorporating strong-field ionization from the excited state. 

PACS numbers: 33.80.Rv,33.20.Xx,33.60.+q,31.15.A- 


I. INTRODUCTION 

The hydrogen molecular ion H J is one of the few systems 
for which the interaction with intense, short laser pulses can 
be simulated truly ab initio, i.e., based on the solution of the 
time-dependent Schrodinger (TDSE) equation Qli without 
further approximations such as, e.g., Born-Oppenheimer or 
Ehrenfest dynamics. Only rotations are usually neglected, 
which is justified for short laser pulses. Despite the simplic¬ 
ity of HJ, its joint electron spectra (JES) for electrons and 
nuclei are intriguingly complex In fact, on top of the 

already complex features in photoelectron spectra from atoms 
ill] there is a nuclear degree of freedom added in (or its 
isotopic sisters). Hence, for any feature observed in a strong- 
field JES at least one question arises: are there vibronic exci¬ 
tations involved? 

Experimental photoelectron ^ectra for and JES for H 2 

have been reported in Refs. simulated ones in Refs. i- 

B). In the multiphoton regime, energy sharing according to 
Eq -h nfiuj = Ee Ep is observed. Here, Eq is the initial 
energy, nfiuj the absorbed photon energy, and Eq, Ep the en¬ 
ergy of the emitted electron and the nuclear kinetic energy 
release (KER), respectively. As Ee + Ep = const., this corre¬ 
lated energy sharing leads to diagonal, straight-line features 
in the E’e, -plane of the JES. At longer wavelengths the 
JES are less simple, especially at low electron energy where 
Coulomb effects are very important, as is well known from 
atomic strong-field ionization M- The diagonal, correlated 
features tend to fade while pronounced oscillations in the ion¬ 
ization probability as function of the electron energy emerge. 
However, also oscillations of the probability for dissociative 
ionization (DI) as function of the KER are observed, which 
have been shown to depend on the initial vibrational state ii. 
One might be tempted to always attribute such variations in 
the DI probability to vibrational excitations. In fact, an inter¬ 
esting application of DI is Coulomb explosion imaging O 
where one strives for reconstructing the initial configuration 
of the nuclei from the KER spectrum after rapid ionization by 
a strong laser field. In this way, e.g., interference structures in 
the KER spectra due to a two-surface population dynamics in 
HJ were observed experimentally We will discuss 

in this paper another mechanism that introduces a modulation 
in the KER. It is based on the oscillatory behavior of the ion¬ 


ization probability as function of the internuclear distance and 
the few-cycle laser pulse duration. 

The paper is organized as follows. In Sec.|II]we start with 
the full quantum H J model and introduce the effect we dis¬ 
cuss in the remainder of this work: the “vertical fringes” (VE) 
in JES, indicating strong variations of the DI yield as func¬ 
tion of the KER but almost independent of the electron en¬ 
ergy. The subsequent sections serve to prove that the VE effect 
is not due to vibrational excitations (Sec. nnii, not due to the 
one-dimensionality of our model (Sec.HVl). and not due to the 
two-center nature of diatomic molecular potentials (Sec. [V]). 
In Sec.[Vll a two-level model, combined with the strong-field 
approximation, is introduced that is capable of qualitatively 
reproducing the VE effect. We conclude in Sec. I VIII and give 
ah the numerical details in the Appendix, in particular on the 
t-SUREE approach for calculating the JES in the various ge¬ 
ometries. 

Atomic units h = tUq = \e\ = 47reo = 1 are used unless 
otherwise indicated. 

II. FULL QUANTUM H+ MODEL 

The Hamiltonian for the two-dimensional H J model reads 

H = + + ( 1 ) 

Both the electronic degrees and the nuclear degrees of free¬ 
dom, i.e., electronic coordinate x G {— 00 , 00 ) and internu¬ 
clear distance i? G (0, cxo), are restricted to the laser polariza¬ 
tion axis, M = 1836 is the proton mass, /i = 2M/{2M + 1), 
and f3 = (M + 1)/M. As we are dealing with a homonuclear 
diatomic molecule the laser field only couples to the electronic 
degree of freedom (a. The velocity-gauge coupling to the 
laser field of vector potential A{t) (in dipole approximation) 
was chosen, with the purely time-dependent A^-term trans¬ 
formed away. 

Eor the interaction between electron and protons we choose 


\/{x - R/2y + e a/(x + i?/2)2 + e 

with the smoothing parameter e = 1. The JES for DI is calcu¬ 
lated via the time-dependent surface flux method (t-SUREE) 
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FIG. 1. (Color online) Computational grid for the full quantum HJ 
model. Coordinates x and R are electronic and internuclear distance, 
respectively. In the blue areas around x — Rj 2 and x — —Rj^ 
the electron is close to one of the protons. The relevant t-SURFF 
boundary for DI (with electrons escaping in positive x direction) is 
given by X = Xdi. The red area indicates the region in which a mask 
function absorbs probability density. 


O (see Appendix [B] for details). Figure [T] depicts the ge¬ 
ometry of the system. Upon time-propagation, probability 
density will pass the surfaces defined by a sufficiently large 
\x\ = Adi, “recorded” there for the calculation of the JES us¬ 
ing t-SURFF, and be absorbed by a mask function thereafter. 
For the t-SURFF approximation we have to assume V^q- ^ 0 
which makes sense if Adi R/2. 

Figure [2] shows correlated spectra for the process of DI in 
few-cycle laser pulses of vector potential 

A{t) = 74osin^ (^) 

for 0 < t < Tp = 27rnc/cc;, calculated by the absolute square 
of (IB6I) . Tp is the pulse duration and Uc the number of laser 
cycles. The initial wave function for the TDSE simulations 
was always the ground state of the Hamiltonian O, which 
has the energy Eq = —0.78. The laser parameters are given 
in the figure caption. 

The main features in the spectra in Fig.[2]are nearly vertical 
and horizontal fringe patterns. Diagonal features indicating 
energy sharing according io Eq ^ nfiuj = Ae + Ep are not 
observed for the laser parameters and the direction of escap¬ 
ing electrons chosen. The modulation of the yield for fixed 
nuclear kinetic energy Ep as function of the electronic kinetic 
energy Ee is well known from laser atom interaction. It can be 
attributed to the interference of several electron paths with dif¬ 
ferent ionization times that lead to the same final electron mo¬ 
mentum (see, e.g., QEIil). Moreover, “direct electrons” 
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FIG. 2. (Color online) Upper panel: JES for DI of HJ in a sin^- 
shaped laser pulse with parameters: ric = 3, A = 800 nm, /peak = 
2.0 • 10^^ W/em^, Up = 0.44. The horizontal line marks the eut- 
off for direet ionization (2Up). The upper JES was ealculated at the 
upper t-SURFF boundary +Adi. Lower panel: ric = 4, /peak = 
1.0 • 10^^ W/em^. The number of VF per nuelear kinetie energy 
interval increases with ric. The lower JES was calculated at the lower 
t-SURFF boundary —Adi. 


and rescattered electrons can be clearly distinguished. The 
“simple man’s” cut-off 2Up is indicated in both panels by a 
horizontal line. The yield due to direct electrons stretches well 
beyond 2Up before it drops down to the level of rescattered 
electrons (approximately four orders of magnitude smaller, 
visible for 0.26 < Ep < 0.47). 

The objective of this paper is to reveal the origin of the mod¬ 
ulation of the yield as function of the nuclear kinetic energy 
Ep. In other words, why is the DI yield strongly suppressed 
for certain proton energies? And why is this suppression al¬ 
most independent of the electronic energy (i.e., why are the 
corresponding fringe patterns almost vertical in Fig. ^ but 
dependent on the pulse duration? Similar modulations have 
been reported in Ref. (I] for simulations starting from a vibra- 
tionally excited H J molecule. The number of VF was found 
to increase with increasing vibrational quantum number u of 
the initial state. However, in our simulations we started from 
z/ = 0 so that the vertical pattern in Fig.[2]does not just refiect 
the probability density of the initial vibrational wave packet. 

On one hand, experience shows that commonly all spectral 
features in strong field ionization can be explained in terms 
of interfering quantum trajectories. On the other hand, the in¬ 
terference of the usual long and short trajectories starting at 
the two nuclear sites (Hi (including potential rescatterin g an d 
the generation of a double-slit type interference pattern 111^ ). 
should depend not only on the internuclear distance but also 
on the electronic energy Eq. Hence the strong suppression 
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FIG. 3. (Color online) Top panel: electronic spectra (each for fixed 
R) for HJ (same laser parameters as in Fig. [51 upper panel). Inlay in 
top panel: spectrum of electrons escaping in the opposite direction. 
Middle panel: occupation of the first excited state, 1— occupation of 
the ground state, and the total ionization probability. Bottom panel: 
energies of ground and first excited state, and their transition dipole 
moment. 

of the DI yield for certain values of Ep but almost indepen¬ 
dent of Ee cannot be explained by such interfering quantum 
trajectories. 

III. FIXED INTERNUCLEAR DISTANCE 

To rule out vibrations as the origin of the VF in Fig.|2]results 
for H J with fixed internuclear distances are discussed now. 

In order to calculate electronic spectra for the H J model 
with fixed inter nuclear distance R the electronic TDSE 
idt'^{x,t) = for the Hamiltonian 

H{t) = -^dl- iA{t)d, + V{x) (4) 

and binding potential 

V{x) = - , ^ - , ^ (5) 

a/(x - Rj2Y + e \/{xE -h e 

was solved for many R in the range where XjR'^ Ep covers 
the relevant KER Ep in Eig.O 

In the top panel of Eig. [3] all these electronic spectra are 
collected for comparison with Eig. [51 upper panel. The over¬ 
all trend is an increasing ionization yield with increasing R 
because of the decreasing ionization potential Ip = |E^o| (see 
bottom panel). The most important insight gained from these 
fixed-i? simulations is that pronounced suppressions of the 


ionization yield are observed for certain internuclear distances 
R l/E’p as well. This proves that vibrational excitation 
cannot be the origin of the VE visible in both Eig.[3]and Eig. [51 

The bottom panel of Eig. [3] shows the energies of the two 
lowest bound electronic states in the H J potential (O vs the 
ionization potential Ip = |E^o|- Fof large internuclear dis¬ 
tances R the two levels are almost degenerate with the ground 
state energy rising asymtotically towards t he grou nd state en¬ 
ergy value for the potential V{x) = —1/ + e. Because of 

this asymptotic degeneracy and the related diverging transi¬ 
tion dipole moment d^i these two states were dubbed “charge 
resonance states”. However, we should stress that the VE in 
the DI yield as discussed in this work occur at smaller dis¬ 
tances than “charge resonance enhanced ionization” (CREI) 

S. 

The middle panel of Eig. [3] shows the total ionization prob¬ 
ability 

P^, = l- E Ax\^{x,T)W ( 6 ) 

J-Xi 

and the occupations of the ground and first excited states after 
the interaction with the laser pulse. The modulations in the 
total ionization probability is less pronounced than in the en¬ 
ergy resolved spectrum, which can be explained by the “left- 
right asymmetry” 0] of the spectra for electrons escaping in 
polarization direction (top panel) and opposite to it (inset in 
top panel). Moreover, the modulations in the energy resolved 
spectrum are not strictly independent of the electronic energy 
Ee, i.e., not perfectly vertical but slightly tilted. 

The fact that ionization probability and bound state occupa¬ 
tions at the end of the pulse oscillate similarly as function of Ip 
{or l/R Ep) suggests that the electronically excited state 
plays an important role in the DI process in few-cycle laser 
pulses. However, the minima in the ionization yield do not 
perfectly coincide with the minima in the excited-state popu¬ 
lation. Eor ionization probabilities smaller than 10“^ the oc¬ 
cupation of the excited state rather oscillates with twice the 
frequency of Pion as function of Ip. Hence, the ionization step 
introduces an additional, nontrivial -dependence. Note that 
the excited state energy |Pi | 0.8 varies little with R so that 

in any case at least about 12 photons are required for ioniza¬ 
tion. Additionally, 2 up to 10 photons are needed to couple 
the initial electronic ground state of energy Eq with the ex¬ 
cited state of energy Pi. Below, in Sec. we reproduce the 
VE qualitatively, using a two-level approximation in combi¬ 
nation with the strong-field approximation (SEA). 


IV. 3D H+ WITH FIXED INTERNUCLEAR DISTANCE 

In order to show that the modulation in the DI yield is not 
an artifact of the low dimensionality of our models the 3D 
molecular ion with fixed ions and aligned along the laser po- 

























4 



l/R (Hartree) 


FIG. 4. (Color online) Same as Fig. [3] but for 3D HJ aligned in 
polarization direction of the laser. 


larization axis was considered. The Hamiltonian 
H = \{-iV + A{t)e,f 


- R/2Y + p2 + i?/2)2 + p2 

is cylindrically symmetric so that the natural choice for the t- 
SURFF boundary is the surface of a cylinder with radius R\ 
and height 2Z\ (see Appendix iDl for details). 

Figure [4] shows spectra for various internuclear distance R 
and electrons escaping in polarization direction. The mod¬ 
ulation of the ionization yield as function of 1/i? is clearly 
visible although at low electron energies the fringes are more 
tilted than in the ID results, making the suppression of the 
yield for certain internuclear distances less electron energy- 
independent. The fringe pattern for the rescattered electrons 
instead is as vertical as in the ID results. Revealing the origin 
of this difference between ID and 3D results requires further 
systematic investigations. In this work we are content with the 
fact that the modulation in the ionization yield as function of 
l/R exists in 3D as well. 


V. SINGLE-CENTER POTENTIAL 

Next we show that the two-center nature of the binding po¬ 
tential is not essential for the observed modulations of the ion¬ 
ization yield, while the existence of an excited state is. To that 
end we consider a Poschl-Teller potential of the form 

V{x) = -- (g) 

8a^ cosh^ [x/{2a)] 

for which the finite number of energy levels of energy = 
— {h — n — 1)^/(8a^), n = 0,1, 2,... < 6 — 1 can be adjusted 
via the parameters a > 0 and 6 > 1. 

First, we aim at mimicking the behavior of ground and ex¬ 
cited state in the molecular model, i.e., Ei = —0.85 is kept 
constant, and Eq covers the range —1.35 < Eq < —1. Figure 
[5] shows the electron spectra collected such that they can be 
directly compared to Fig. O The VF are there, proving that 
they are not due to a two-center interference. 
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FIG. 5. (Color online) Same as Fig.[3]but for a Poschl-Teller potential 
(|8]) with ground and excited state energy tuned close to the molecular 
case. The laser intensity was increased to /peak = 3.0 • 10^^ W/cm^ 
in order to have an ionization yield similar to the molecular models. 
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FIG. 6. (Color online) Same as Fig.[5]but for a Poschl-Teller potential 
(El) supporting only a single bound state. The intensity was increased 
further to /peak = 6.0-10^^ W/cm^. 


Second, Fig. [6] shows the case of a Poschl-Teller potential 
with a single bound state only. The intensity was increased 
to /peak = 6.0 • 10^^ W/cm2 to Compensate for the decreasing 
ionization in the narrower and deeper potential. Each individ¬ 
ual photoelectron spectrum looks standard “SFA-like”. Both 
the VF and oscillations in the occupation of the groundstate 
at the end of the laser pulse are absent. This substantiates our 
assertion that the occupation of an excited state is crucial for 
the modulation of the (dissociative) ionization yield. 


VI. TWO LEVELS + SEA 

As long as the ionization probability is small we may 
model the occupation of the ground and first excited states 
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FIG. 7. (Color online) Upper panel: photoelectron spectra for the 
model of Fig. [3] calculated using the two-level SFA. Lower panel: 
occupation of the first excited state as calculated from the TDSE, the 
numerical solution of Eqs. (fTol) and (ITTI) . and the analytical result 

d. 


by a simple two-level model. Plugging the ansatz \'ip{t)) = 
ao(t)IV^o) + into the TDSE in length gauge 

idt'ip{x,t) = +E{t)x+ V{x)^'tp{x,t) (9) 


one finds the well-known equations of motion for the density 
matrix elements pij = a*aj, i = 0,1, j = 0,1, 

Poo = idoiE{t){pio — pIq) (10) 

pio = idoiE{t){poo — pii) iAE'pio (11) 

where AE = Ei — Eq, doi = { 2 po\x\ 2 pi) (assumed real), 
Poi = Pio^ Pii = 1 — Poo- As we are interested in few-cycle 
pulses and the transient dynamics induced by them we can¬ 
not apply the rotating wave approximation, and a dressed or 
Floquet state approach does not make sense either. Instead, in 
the bottom panel of Fig. [7] the density matrix element pn at 
the end of the laser pulse t = from the numerical solution 
of the two-level model Eqs. ([TOb and ([TT1) [initial conditions 
poo{t = 0) = 1 and pii{t = 0) = 0] is compared to the occu¬ 
pation of the first excited state from the numerical solution of 
the full TDSE. The agreement is very good apart from a shift 
along the ionization potential axis. This shift is caused by ne¬ 
glecting higher excited states and the coupling to the contin¬ 
uum in the two-level model. We checked that for lower field 
strengths (where even less ionization occurs and other excited 
states are negligibly populated) the agreement improves. 

The results of the two-level system can be used to model 
ionization as well. In “standard” SFA only a single bound 
state (besides the continuum states of momentum k) is con¬ 
sidered, and depletion of its population is neglected (see, e.g.. 
Ref. 13]). Instead, we plug the modulus of the occupation 
= \/pii(t) of the first excited state into the SFA am¬ 
plitude for direct ionization. 


\cil, sfa(^,^)|^ 




( 12 ) 


Here, dik{k) = {'ipi\x\k), and we neglect the amplitude for ionization from the ground state because its contribtion is several 
orders of magnitude smaller than the contribution of the first excited state. This is fortunate, as otherwise the phases of the 
complex ao{t) and ai{t) matter and should be calculated from a full SFA with two bound states and all relevant bound-bound 
and bound-continuum couplings. The top panel in Fig. [7] shows the collected electronic spectra for the molecular potential with 
fixed protons of Sec. Hill calculated using our simplified two-level SFA. The dipole moments and the eigenenergies Eq and 
El were calculated numerically from the TDSE data of Sec. [III| (see the bottom panel of Fig. 0 for the i7-dependent energies 
and transition dipole moment d^i). Figure [7] shows that the simple two-level SFA reproduces the -dependent features in the 
ionization probability qualitatively. In particular, the correlation between the oscillations of the excited-state population at the 
end of the pulse as function of Ip and the oscillations in the ionization probability with only half the frequency is as observed in 
the TDSE results of Sec.HIIl 

Employing pu < poo = 1 - Pii — 1 in (HI]), Eqs. dlOl) and (II 111 can be solved analytically, leading to 




— r dr'£;(r') 

w J 


T = n = AE/uj. 


(13) 


Using E{t) = —dtA{t) with the vector potential ([3]) yields the occupation of the first excited state at the end of the sin^ pulse 

AQdoiUJ^{n‘^uj‘^ — -h 3AE‘^n‘^)AE sin 
= -?-u-77-^—77-—7 

00-00-1-~ ~ [^cOO- -j- 00) [ricOO-i. -h 00) 


where co’+ = cc + AE and oo-^ = 00 — AE. Expanding this expression in the small parameter 7 ^ = oojAE = 1/n gives 

dA^d^i sm^{7rAEnc/oo)p^ 


( 15 ) 
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The occupation of the first excited state after the pulse thus 
decreases as increases. Hence the observed VF in the (dis¬ 
sociative) ionization yield are a few-cycle effect. Moreover, 
inspection of the sine’s argument in (O shows that the fre¬ 
quency of the oscillation depends on the number of cycles ric. 
The higher ric the more oscillations within a given Ain¬ 
terval. This is in agreement with the TDSE results in Fig. [2| 
where the 4-cycle laser pulse was found to generate more VF 
than the 3-cycle pulse. As Ei is almost constant in the 
model the oscillations in Fig. [ 7 ] are of almost constant period 
when plotted vs Ip = jE'o I • 

Note that pn is very sensitive to the pulse shape. In fact, 
for a Gaussian pulse pn according (fTSl) —with the integra¬ 
tion limits stretched to ±oo —becomes the (modulus squared) 
Fourier transform of a Gaussian, which is a Gaussian and thus 
does not oscillate with AE" and ric (full-width half maximum) 
at all. 


VII. CONCLUSION 

Numerical simulations of the dissociative ionization pro¬ 
cess in for short laser pulses reveal patterns of vertical 
fringes in the joint energy spectra, i.e., strong variations of the 
yield as function of the ion energy that are almost independent 
of the electron energy. Identifying the kinetic energy release 
with the inverse internuclear distance, the effect is also found 
in calculations with fixed ions, ruling out vibrational excita¬ 
tions as its origin. Instead, ionization proceeds via the first 
excited electronic state. In few-cycle pulses the population of 
the first excited state depends strongly on the number of cy¬ 
cles and the pulse shape in general. The vertical fringes in the 
continuous dissociative ionization spectra are clearly corre¬ 
lated with the population of the first excited state at the end of 


the pulse, as qualitatively reproduced using a simple two-state 
model combined with the strong-field approximation. 

The observed effect relies on the ultrashort, transient dy¬ 
namics in few-cycle laser pulses and not on resonances, spe¬ 
cially chosen detunings, or interference. In fact, in the limit 
of long laser pulses the vertical fringes disappear and one 
approaches—depending on the laser frequency—either or¬ 
dinary non-resonant multiphoton or tunneling ionization, or 
well-known resonance-enhanced multiphoton ionization. 
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Appendix A: Numerical details 

The TDSE was solved numerically by propagating the 
wavefunction with the Crank-Nicolson time propagator. The 
wavefunction and the potentials were discretized on a Carte¬ 
sian grid with the spatial derivatives in the Hamiltonian ap¬ 
proximated by finite differences. An iterative block Gauss- 
Seidel method and the Thomas algorithm were applied for 
the solution of the linear system of equations of the Crank- 
Nicolson method in the two and one dimensional case, re¬ 
spectively. The initial ground-state wavefunctions for the time 
propagation and the first excited state were obtained by the 
shift-invert method EH. For the numerical solution of the 
cylindrically symmetric Hamiltonian (| 7 ]) the coordinate trans¬ 
formation ^ was used Ei. Numerical parameters for 

the TDSE simulations are summarized in Table HI 


Appendix B: t-SURFF for the HJ model 


Assuming that V^q- in (O can be neglected for x > Xdi the wavefunction there separates in the form '0/c(x, t)(j)p{R) e 
where (j)p{R) are the solutions of the Coulomb scattering problem 

(tn + 2) 0p(i?) = , Tn = -^dl, (Bl) 

and 

1 

'0/c(x, t) = (27r)“^/^ a{t) = - / dt'[k‘^ E 2kA{t')] (B2) 

2 Jo 


are Volkov wavefunctions. 

The DI amplitude (restricted to the electrons escaping in positive direction) is approximated by the integral 

{'^{T,)\DI)c^aMk,p) = {^{T)\e{x-Xj,,)\p{T))\k{T)) = JdR J dx R,T)Mx,T)MR) ■ (B3) 

x>Xy)i 

This expression is not yet useful for practical purposes because T needs to be large enough to allow the slow electrons arriving 
in the region x > Xdi. On the other hand, the fast electrons need to be kept on the grid as well, necessitating a huge grid size. In 
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order to avoid large grids the t-SURFF method [H, [TsI] was adapted to the problem at hand. Writing the right hand side of (IB3I) 
as a time integral we obtain 


aDi { k , p ) = (^(O)|0(x - Xdi)|p(0))|/c(0)) + [ dt dt {^{ t )\ e{x - Xj^i) \ p { t ))\ k { t )) . (B4) 

Jo 

For sufficiently large Xdi and a bound initial state |^) we have (^(0) |0(x — Xdi) |p(0)) |/i:(0)) 0. Employing the TDSE with 

0 yields 


CiDl { k , p ) aDI,t-SURFF(^,p) = M (^(01 

Jo 


- ^Di) \pit))\Ht)) ■ (B5) 


Only terms of the Hamiltonian containing derivatives with respect to x contribute in the commutator, leading to 

pT ^oo 
/O ^0 


pT poo 

ciDi,t-suRFF{k,p) = dt dR PA{t)^''{X]:,i,R,t)'ipk{^Di,t) 
Jo Jo ^ 


- -^{'^'^{Xm,RR)dx'ipk{xR)\x=x^i - 'ipk{XmR)dx'^''{x,RR)\x=x^,) (l)p{RR). (B6) 

The scattering states (l)p{R) were used as implemented in the GNU Scientific Library (GSL) ll^ . 


Appendix C: t-SURFF for ID calculations 

For the one dimensional systems the probability amplitude for ionization with final electron momentum k is approximated as 

(Cl) 

Again, only electrons escaping in positive x direction, passing the t-SURFF boundary Xi, are considered, and the binding 
potential is neglected for distances x > Xi. 


Appendix D: t-SURFF for cylindrically symmetric system 


The probability amplitude for an electron escaping with a momentum k = kpSx + kzCz can be approximated by the integral 


(Dl) 


{k{T)\^{T)) ^ dt JdVdt(^Mp,z,ty{o{Ri-p)o{z-Zi) + e{Ri-p)e{-z-Zi) + e{p-Ri))^{p,z,t)'). 

The integral (IDlb can be divided into three terms which are evaluated separately. Using the TDSE, the first term reads (dropping 
the arguments of 'ipk and 

si{kp, ky = dt I dVdt (ymRi - p)0{z - Z,)^'^ 

= -^ £dt I dv^k [p~"dppdp + dl + 2iA{y)d,, 0(i?i - p)e{z - z^)] ^ 

= I dtj dvrk {o{z - Zi) {^dp5{Ri -p) + 2{dp^)5{Ri - p) + p-H5{Ri - p)) 

-e{Ri - p) yid.Siz - Zi) + 2{d,-^)5{z - Zi) + 2iA{t)-^5{z - Zi)) ) 

• pT ^2 77 / pOO 

^ ^ Jo Jo ^xJz + 2pV’fc(5p5') + 'Ipk'^) 

- £' pdp {-dyrk^) + 2rk{d.^) + 2iA{t)rk^) . 


\p=Ri 


(D2) 
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The second term is, analogously, 

rT r 27 T / n-Zi 


S2{kp,k:,) = ^ dt dz{-dp{p'tpk'^) + 2pi}l.{dp'^)+ipl.'^)\p=R, 

+ r pdp{-d,{rk^) + 2rk{d,^) + ‘liAitWk^) U=_z, 

Jo 


, (D3) 


and 


ssikp, kO = dtjdVdt[rkO{p - 

= ^ dt dp dz{-dp{p^^k)+ ‘^P^ki9p^)+^k^)\p=R- 
^ Jo Jo j-oo 

Inserting the Volkov wavefunction 

^fc(p, t) = (27r)-3/2 e-i»(t)+ik.z+ik,pcos^^ ^ 1 / [^2 2k^A{t’)] 

2 Jo 

and collecting all integrals over 2 ; in S 3 such that s'^ + S 2 + S 3 = si + S 2 + S 3 yields 

rT nRi 




s'i{kp,kz) = --{27r) e 


[ ' pdpJo{kpp){ik,^ + 2iA^) 

Jo Jo 


z=Zi 5 


S2{kp,k0 = i/2g*fczZi ^ ^^gia(t) / pdpJo{kpp) {ikz^ + + 2iA'^)\z^_Zi, 

^ Jo Jo 

s'^{kp,kz) = ^(27r)-i/2i?i [ dte^^d) f ' dze-^’^^^ (^kpMkpRj) + JoikpRi)dp^) \p=R, 

^ Jo J- 


and 


Jo J-Zi 

where Jq and Ji are Bessel functions of the first kind. Denoting the Fourier transform 

rZi 

^{p, kz,t) 


l-Zi 


dze-^’‘^^^{p,z,t) 


and the Hankel transform 

rR\ 

'^{kp,z,t)= dppJo{kpp)'^{p,z,t) 
Jo 

the approximation of the probability amplitude reads 


(D4) 

(D5) 

(D6) 

(D7) 

(D8) 

(D9) 

(DIO) 


Si + S 2 + S 3 = 


2V^ 


fT 

/ dt {ikz + dz + 2iA) ^{kp, z, t)U=z, - {ik^ + dz + 2iA) ^{kp, z, t)\z=-z, 

Jo 

+ Ri{kpJi{kpRi) + Jo{kpRi)dp)'^{p,kz,t)\p=Ri) ■ (Dll) 


The GSL Q was used for the Hankel transform and the Bessel functions. 

In order to suppress spurious effects introduced by the finite time T in the t-SURFF time integrals a Hanning window 


H(t) = 

was multiplied to the integrands (1^ - (IcTTi . and (ITTrni 


1 iff < T/2 

[1 — cos(27rf/r)]/2 lit > T/2 


(D12) 
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